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ABSTRACT 

The core rotation rates of massive stars have a substantial impact on the nature of core-collapse 
supernovae and their compact remnants. We demonstrate that internal gravity waves (IGW), excited 
via envelope convection during a red supergiant phase or during vigorous late time burning phases, can 
have a significant impact on the rotation rate of the pre-SN core. In typical (10 M 0 < M < 20 M 0 ) 
supernova progenitors, IGW may substantially spin down the core, leading to iron core rotation periods 
-Pmin.Fe > 30 s. Angular momentum (AM) conservation during the supernova would entail minimum 
NS rotation periods of P m i n ,NS ^ 3 ms. In most cases, the combined effects of magnetic torques and 
IGW AM transport likely lead to substantially longer rotation periods. However, the stochastic influx 
of AM delivered by IGW during shell burning phases inevitably spin up a slowly rotating stellar core, 
leading to a maximum possible core rotation period. We estimate maximum iron core rotation periods 
Of P m ax,Fe < 5 X 10 3 S in typical core-collapse supernova progenitors, and a corresponding spin period 
of Pmax.NS 1$ 500 ms for newborn neutron stars. This is comparable to the typical birth spin periods 
of most radio pulsars. Stochastic spin-up via IGW during shell O/Si burning may thus determine the 
initial rotation rate of most neutron stars. For a given progenitor, this theory predicts a Maxwellian 
distribution in pre-collapse core rotation frequency that is uncorrelated with the spin of the overlying 
envelope. 


1. INTRODUCTION 

Rotation is a key player in the drama that unfolds upon 
the death of a massive star. The angular momentum 
(AM) contained in the iron core and overlying layers de¬ 
termines the rotation rate at core collapse (CC), which 
could have a strong impact on the dynamics of CC and 
the subsequent supernova (see e.g MacFadyen & Woosley 
1999; Woosley & Heger 2002, 2006; Yoon ct al. 2006). 
Rotation may help determine the nature of the compact 
remnant, which could range from a slowly rotating neu¬ 
tron star (NS) to a millisecond magnetar or rapidly ro¬ 
tating black hole (see e.g. Heger et al. 2000; Heger et al. 
2005). The former may evolve into an ordinary pulsar, 
while the latter two outcomes offer exciting prospects 
for the production of long gamma-ray bursts (GRB) and 
superluminous supernova. In each of these phenomena, 
rotating central engines are suspected to be the primary 
source of power (Woosley 1993; Kasen & Bildsten 2010; 
Metzger et al. 2011). 

Despite rotation being recognized as an important 
parameter controlling the evolution of massive stars 
(Maeder & Mcynet 2000), little is known about the rota¬ 
tion rates of the inner cores of massive stars nearing CC. 
The best observational constraints stem from measure¬ 
ments of the rotation rates of the compact remnants fol¬ 
lowing CC. For instance, a few low-mass black hole X-ray 
binary systems have been measured to have large spins 


that can only be accounted for by high spins at birth (Ax- 
elsson et al. 2011; Miller et al. 2011; Wong et al. 2012). 
However, the rotation rates of young NSs show little evi¬ 
dence for rapid rotation (P < 10 ms) at birth. The most 
rapidly rotating young pulsars include PSR J0537-6910 
(P = 16 ms) and the Crab pulsar (P = 33 ms), whose 
birth periods have been estimated to be Pi < 10 ms (Mar¬ 
shall et al. 1998) and Pi ~ 19 ms (Kaspi & Hclfand 2002), 
respectively. Many young NSs appear to rotate much 
more slowly, with typical periods of hundreds of ms (l.ai 
1996; Gotthelf et al. 2013; Kramer ct al. 2003; Dingel 
et al. 2015). In general, pulsar observations seem to in¬ 
dicate a broad range of initial birth periods in the vicinity 
of tens to hundreds of milliseconds (Faucher-Giguere & 
Kaspi 2006; Popov et al. 2010; Popov & Turolla 2012; 
Gullon et al. 2014). Hence, rapidly rotating young NSs 
appear to be the exception rather than the rule. 

Theoretical efforts have struggled to produce slow rota¬ 
tion rates. In the absence of strong AM transport mech¬ 
anisms within the massive star progenitor, NSs would 
invariably be born rotating near break-up (Heger et al. 
2000). Heger et al. (2005) and Suijs ct al. (2008) ex¬ 
amined the effect of magnetic torques generated via the 
Tayler-Spruit (TS) dynamo (Spruit 2002), and found 
typical NS spin periods at birth (assuming AM conserva¬ 
tion during CC and the ensuing supernova) of P ~ 10 ms. 
Wheeler et al. (2014) implemented magnetic torques due 
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to MRI and the TS dynamo, and were able to reach iron 
core rotation rates of P c ~ 500 s, corresponding to NS 
spin periods of P ~ 25 ms. These efforts are promis¬ 
ing, but the operation of both mechanisms within stars 
has been debated (e.g. Zahn et al. 2007), and theoretical 
uncertainties abound. 

Recent asteroseismic advances have allowed for the 
measurement of core rotation rates in low-mass red gi¬ 
ant stars (Beck et al. 2012, 2014; Mosser et al. 2012; 
Deheuvels et al. 2012, 2014). In these stars, the core ro¬ 
tates much faster than the surface and one cannot assume 
nearly rigid rotation as suggested in jpruit & Phinney 
(1998). However, the cores of low-mass red giants rotate 
much slower than can be explained via hydrodynamic 
AM transport mechanisms or magnetic torques via the 
TS dynamo (Cantiello et al. 2014). If similar AM trans¬ 
port mechanisms operate in more massive objects, this 
suggests that the pre-collapse cores of massive stars may 
rotate slower than predicted by many previous theoreti¬ 
cal investigations. 

Internal gravity waves (IGW) constitute a powerful en¬ 
ergy and AM transport mechanism in stellar interiors. 
Several studies (Kumar & Quataert 1997; Zahn et al. 
1997; Kumar et al. 1999; Talon et al. 2002; Talon & 
Charbonnel 2003, 2005, 2008; Charbonnel & Talon 2005; 
Denissenkov et al. 2008; Fuller et al. 2014, hereafter F14) 
have found that convectively generated IGW can redis¬ 
tribute large quantities of AM within low-mass stars. 
IGW may partially account for the rigid rotation of the 
Sun’s radiative interior and the slow rotation of red giant 
cores, although magnetic torques are also likely to be im¬ 
portant (Denissenkov et al. 2008, F14). IGW may also 
be important in massive main sequence stars. Rogers 
et al. (2012) and Rogers et al. (2013) have found that 
convectively generated IGW can alter the spin-rate of 
the stellar photosphere, while Lee et al. (2014) propose 
that convectively generated IGW in B[e]-type stars may 
instigate outbursts that expel mass into the decretion 
disk. 

The AM redistribution arising from convectively ex¬ 
cited IGW stems primarily from two physical effects. 
The first, more commonly studied effect, is wave filtering 
via differential rotation, which stems from an IGW-mean 
flow interaction. The second, less commonly studied ef¬ 
fect, is stochastic variations in the AM flux carried by 
IGW, which arises because the IGW are excited by the 
stochastic motions of turbulent convection. We shall see 
that both mechanisms are important at different phases 
of the lives of massive stars. 

Convectively excited IGW have a particularly strong 
influence on the evolution of massive stars nearing CC. 
Indeed, after core carbon exhaustion, waves are the most 
effective energy transport mechanism within radiative 
zones, as photons are essentially frozen in and neutri¬ 
nos freely stream out. In two recent papers, Quataert 
& Shiode (2012) and Shiode & Quataert (2014) (here¬ 
after QS12 and SQ14) showed that the prodigious power 
carried by convectively excited waves (on the order of 
10 10 L e during Si burning) can sometimes unbind a large 
amount of mass near the stellar surface, and may sub¬ 
stantially alter the pre-collapse stellar structure. IGW 
are ubiquitous in simulations of late burning stages 
(Meakin & Arnett 2006, 2007; Meakin & Arnett 2007), 


although existing simulations have not quantified their 
long-term impact. 

In this paper, we examine AM transport due to con¬ 
vectively excited IGW within massive stars, focusing pri¬ 
marily on AM transport during late burning stages (He, 
C, O, and Si burning). We find that IGW are generally 
capable of redistributing large amounts of AM before CC 
despite the short stellar evolution time scales. In the 
limit of very efficient prior core spin-down, we show that 
a stochastic influx of AM via IGW sets a minimum core 
rotation rate which is comparable to the broad distribu¬ 
tion of low rotation rates (P < 500 ms) observed for most 
young NSs. We also examine whether IGW can prevent 
the core from spinning up as it contracts during stellar 
evolution. We find that IGW may spin down the core 
substantially, although spin-down via magnetic torques 
is likely required to reproduce the observed pulsar pop¬ 
ulation. 

Our paper is organized as follows. In Section 2, we 
describe our massive star models, the generation of IGW 
during various stages of stellar evolution, and the AM 
they transport. In Section 3, we consider whether IGW 
can stochastically spin up a very slowly rotating core, 
attempting to determine a maximum core rotation pe¬ 
riod. In Section 4, we investigate whether the IGW can 
spin down the cores of massive stars, attempting to de¬ 
termine a minimum core rotation period. In Section 5, 
we conclude with a discussion of our results and their im¬ 
plications for CC, supernovae, and the birth of compact 
objects. 

2. CONVECTIVELY EXCITED IGW 

IGW are generated by convective zones and propagate 
into neighboring stably stratified regions, carrying en¬ 
ergy and AM. To estimate energy and AM fluxes car¬ 
ried by IGW, we use techniques similar to those of F14, 
QS12, and SQ14. We begin by constructing a sequence 
of stellar models using the MESA stellar evolution code 
(Paxton et al. 2011, 2013). In what follows, we focus on 
a M = 12 Mq, Z = 0.02 model that has been evolved 
to CC, although we also discuss a more massive (35 Mq) 
model in Section 3. Details on the models can be found 
in Appendix A. For our purposes, the most important 
model outputs are the local heat flux, convective rnach 
numbers, and life time of convectively burning zones. As 
in SQ14, we find these quantities correlate most strongly 
with the helium core mass. Stellar models of larger zero- 
age main sequence (ZAMS) mass or with more mixing 
(due to overshoot or rotation) tend to have a higher He 
core mass and may exhibit different wave dynamics than 
our fiducial model. Our main goal here is simply to pro¬ 
vide a rough estimate of IGW AM fluxes for a typical 
low-mass (M < 20 Mq) progenitor of a NS. 

A full understanding of AM transport by IGW should 
include the combined effects of waves emitted from each 
convective zone. For simplicity, we focus on cases in 
which a convective shell overlies the radiative core, irra¬ 
diating it with IGW. These convective shell phases typi¬ 
cally occur after core burning phases and thus have the 
final impact for a given burning phase. We use mixing 
length theory (MLT), as described in F14, to calculate 
IGW frequencies and fluxes. Encouragingly, our MLT 
calculations yield convective velocities and Mach num¬ 
bers similar to those seen in simulations. We find convec- 
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Fig. 1.— Kippenhahn diagram of our 12 Mq stellar model, with 
the x-axis showing the time until CC. Light blue shaded regions are 
convective. Solid colored lines show helium, carbon, oxygen, sili¬ 
con, and iron core masses. We have labeled convective core burning 
phases with the element being burned in each phase. Dashed verti¬ 
cal lines show the locations (from left) of our helium core burning, 
carbon shell burning, oxygen shell burning, and silicon shell burn¬ 
ing models used in Figures 2 and 5. The mass coordinate of the base 
of the convective shell, from which IGW propagate downwards, is 
labeled by a circle on each dashed line. 


tive velocities during Si burning very similar to those of 
Couch et al (2015). Our 12Af 0 O burning model yields 
convective velocities and turnover frequencies a factor of 
~ 2 smaller than those seen in simulations of O burn¬ 
ing shells (e.g., Meakin & Arnett 2006; Meakin & Arnett 
2007; Meakin & Arnett 2007; Arnett ct al. 2008). How¬ 
ever, this is likely due to the larger mass progenitor of 
their model, as we find higher mass models (~ 20 Mq 
as used in the works above) yield O burning properties 
closer to those seen in the simulations. We proceed with 
our MLT results and expect that realistic wave frequen¬ 
cies may differ from those used here by a factor of <2. 

Figure 1 shows a Kippenhahn diagram for our stellar 
model, and Figure 2 shows the density (p), mass [M(r)\, 
and Brunt-Vaisala frequency (TV) profiles of our model 
during important convective shell phases. The first con¬ 
vective shell phase occurs during He-core burning, at 
which point the star has evolved into a red supergiant. At 
this stage, IGW are generated at the base of the surface 
convection zone and propagate toward the He burning 
core. We have also shown profiles during shell C burn¬ 
ing, O burning, and Si burning, when the radiative core 
contains a mass of M c ~ 1 Mq and is being irradiated 
by IGW generated from the overlying convective burn¬ 
ing shell. The basic features of each of these phases is 
quite similar, the main difference is that more advanced 
burning stages are more vigorous but shorter in duration. 
We find that the characteristics of the convective burn¬ 
ing shells (convective luminosities, turnover frequencies, 
mach numbers, and lifetimes) are similar to those listed 
in QS12 and SQ13, although the shell burning phases are 



Fig. 2. — Internal structure of our 12 Mq stellar model, at the 
phases labeled in Figure 1. At all stages shown, the star is a red 
supergiant with radius R 10 3 Rq- The He-burn stage is core 
helium burning, while the C-burn, O-burn, and Si-burn stages are 
shell carbon, oxygen, and silicon burning phases, respectively. Ver¬ 
tical shaded bars are drawn near the base of the surface convection 
zone for the He-burn stage, and near the base of the convective shell 
burning region of the other stages. Top: Density profiles p of the 
inner regions of our models. The envelope density profile is similar 
for each model, but central densities vary by orders of magnitude. 
Middle: Enclosed mass profiles. For C/O/Si burning, the mod¬ 
els are chosen such that the mass internal to the convective shell 
burning region is M c ~ 1.2Mq. Bottom: The Brunt-Vaisala fre¬ 
quency, N, of our models. We have also plotted horizontal dashed 
lines at the convective turnover frequencies uj c below the relevant 
convective zone. IGW propagate in the regions where lj c < N , 
shown by the horizontal extent of the dashed lines. 


generally more vigorous and shorter-lived than the core 
burning examined in SQ13. Table 1 lists some of the 
parameters of our convective zones. 

The total energy flux carried by waves emitted from 
the bottom of the convective zone is of order 

ml c <E< M 5/8 L c (1) 

(Goldreich & Kumar 1990; Kumar ct al. 1999), where M. 
is the convective Mach number (defined as the ratio of 
MLT convective velocity to sound speed, v c /c s ), and L c 
is the luminosity carried by convection near the base of 
the convective zone. Many previous works have used the 
left-hand side of equation 1 as an estimate for the IGW 
energy flux, although Lecoanct & Quataert (2013) argue 
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TABLE 1 


Burning Phase 

r c 

(k 

:m) 

^shell (s) 

t waves (s) 

M 

L 

c ( L @) 

Ldc 

C 

rad/s) 

He Core Burn 

1.6 

X 

10 v 

4 x 10 13 

2 

x 10 s 

0.06 

6 

x 10 4 

3 

X 

10- 

-6 

C Shell Burn 

9.7 

X 

10 3 

3 x 10 s 


10 6 

0.002 

3 

x 10 8 

4 

X 

10- 

-3 

O Shell Burn 

3.6 

X 

10 3 

4 x 10 6 


10 s 

0.004 

8 

x 10 9 

2 

X 

10- 

-2 

Si Shell Burn 

1.7 

X 

10 3 

7 x 10 3 

2 

TO 

o 

X 

0.02 

2 

x 10 12 

4 

X 

10- 

-1 


Properties of convection during the late burning phases shown in Figures and 1 and 2. Here, r c is the radius at the base of the 
convection zone in consideration, T s h e ii is the duration of the burning phase, t wa , vea is a wave spin-up timescale (Equation 4), M 
is the convective Mach number, L c is the luminosity carried by the convective zone, and oj c is the angular convective turnover 
frequency. 


TABLE 2 


Burning Phase 

r c (km) 

^shell ( s ) 

t waves (s) 

M 

Lc (To) 

lj c (rad/s) 

C Shell Burn 

2.0 x 10 4 

7 x 10 s 

10 4 

0.004 

9 x 10 9 

3 x 10 -3 

O Shell Bum 

4.6 x 10 3 

9 x 10 5 

5 x 10 3 

0.008 

2 x 10 11 

4 x 10“ 2 

Si Shell Burn 

2.3 x 10 3 

5 x 10 3 

4 x 10 2 

0.02 

6 x 10 12 

3 x lO” 1 


Same as Table 1, but for our 35 Mq model. We do not examine core He burning for this model because it is not a red supergiant 
and does not have a convective envelope at this stage. 


that a more accurate estimate may be E ~ A 4 5 ^ 8 L C , 
which is larger by a factor of A4 -3 / 8 . We consider the 
left-hand side of equation 1 to be a lower limit for the 
wave flux, and the right-hand side to be an upper limit. 

For shell burning, this energy flux is dominated by 
waves with horizontal wave numbers and angular fre¬ 
quencies near to ~ max(r c /iT c , 1), and ui ~ ui c , respec¬ 
tively. Here, H c and r c are the pressure scale height and 
radial coordinate near the base of the convective zone, 
and we define the angular convective turnover frequency 
as 


w c 


TTV C 

a.W c ’ 


( 2 ) 


where aH c is the mixing length. Our models use a = 1.5. 
For the convective shells we consider, H c ~ r c and we ex¬ 
pect waves of low angular degree to be most efficiently ex¬ 
cited. The characteristic AM flux carried by these waves 
is 


j 



UI 


( 3 ) 


Turbulent convection generates waves with a spectrum 
of azimuthal numbers m and angular frequencies u>. The 
values given above are characteristic values which domi¬ 
nate the AM flux. The waves carrying the most AM flux 
sometimes damp before they reach the core, and might 
not be able to affect the spin of the core. Then the waves 
with to and ui would not dominate the AM flux to the 
core; instead, other waves in the turbulent spectrum be¬ 
come important (see Section 4). 

As a first check to see if IGW can have any affect on 
the spin of the core of the star, we assume all waves 
can propagate to the core. We suppose that IGW could 
be important for the spin evolution if they are able to 
carry an amount of AM comparable to that contained in 
a young NS, which contains J^s ~ 10 48 gcm 2 s _1 for a 
rotation period of Pns = 10 ms. Then the characteristic 
timescale on which waves could affect the AM of the core 
is 


t 


waves 


Jns 

~T 


( 4 ) 


Table 1 lists shell burning lifetimes T s h e n and wave spin- 
alteration time scales f waves, evaluated using E oc A4, 
to = 1 and ui = lo c . In all phases, fwaves <C T s h e ii for 
our model, indicating that waves may be able to have 
a substantial impact on the spin rate of the core. We 
examine this impact in the following sections. 

3. STOCHASTIC SPIN-UP BY INTERNAL GRAVITY 
WAVES 

The core of a massive star undergoing shell burning 
is irradiated by IGW, and so it will generally contain a 
non-zero amount of AM and have a non-zero rotation 
frequency. Our goal here is to calculate the net AM de¬ 
posited by IGW in an initially non-rotating core, and 
therefore to calculate a minimum core spin rate follow¬ 
ing each shell burning phase. 

Our calculation is based on two key ideas. First, we 
show in Appendix B that the core will absorb all of the 
AM carried inward by IGW because the waves become 
non-linear due to geometric focusing near the center of 
the star. This process has been extensively discussed 
in the literature and shown to occur in hydrodynamical 
simulations (Press 1981; Barker & Ogilvie 2010; Barker 
2011). Section 4 also demonstrates that the IGW gener¬ 
ated by late shell burning (C burning and beyond) are es¬ 
sentially unattenuated by radiative diffusion or neutrino 
damping as they propagate inward (see Appendix C for 
additional discussion on the effects of radiative diffusion). 
IGW emitted from a convective shell will therefore prop¬ 
agate toward the center of the star, non-linearly break, 
and deposit their AM within the core. The total AM of 
the core at the end of a shell burning phase will there¬ 
fore be a superposition of the AM deposited by IGW 
throughout the shell burning phase. 

The second key point is that IGW are generated by 
the turbulent, stochastic motions of the convective burn¬ 
ing shell. The AM flux carried by IGW therefore varies 
stochastically on a convective turnover time scale. For 
sufficiently long time scales we expect a nearly zero net 
AM flux because equal amounts of prograde and retro¬ 
grade waves will be generated, but on convective turnover 
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time scales the instantaneous AM flux will be non-zero. 
The net AM flux deposited by IGW during the finite 
duration of the shell burning phase will therefore also 
be non-zero. This discreteness of the AM flux is particu¬ 
larly important in late stage shell burning phases because 
they are short lived compared to main sequence burning 
phases. We note that stochastic IGW AM transport has 
also been invoked in Rogers et al. (2012) and Lee et al. 
(2014) to modify the AM contained in the atmospheres 
of stars on the main sequence. 

A complication is that IGW absorption at critical lay¬ 
ers could occur if IGW generate large shear near the 
center of the star due to its low moment of inertia. 
These critical layers can prevent IGW from entering the 
core if they can migrate outward to the core boundary, 
which occurs on the time scale t s hear = / c w c / J. We find 
Ghear ~ 2 x 10 4 s for Si shell burning, using the conser¬ 
vative estimate from equation 1. Although IGW may 
generate significant shear within the T s h e ii ~ 7 x 10 3 s Si 
shell burning life time, the shear will not encompass the 
entire core and prevent further influx of IGWs. We do 
not attempt to calculate the precise distribution of AM 
within the core (i.e., the core rotation profile, which will 
depend on the non-linear wave breaking dynamics and 
AM redistribution within the core), but only calculate 
the net AM deposited by IGW within the core. 

To quantify our arguments, we consider a non-spinning 
core being irradiated by IGW emitted from an overlying 
convective shell. We define an IGW “packet” as the IGW 
generated over a convective turnover time r c = 27r/w c . 
Since the convective motions are only correlated over 
~t c , each wave packet will have an AM vector uncorre¬ 
lated with the previous wave packet. The spherical sym¬ 
metry of the background structure implies that these AM 
vectors will be randomly oriented in space. From equa¬ 
tion 3, it follows that each wave packet will have an AM 
vector of approximate length 


u> 


2mir 


ML C 


(5) 


The expected magnitudes of the x , y, and z components 
of the AM vector are J w /\/3. Since our goal is to find a 
minimum rotation rate, we use fh = 1 to minimize the 
AM carried by each wave packet. 

As IGW packets deposit their AM in the core, the total 
core AM in each direction exhibits a random walk. Dur¬ 
ing each shell burning phase, the convective shell emits 
roughly TV wave packets, with 

TV ~ W ^ he11 , (6) 

27T 

where T s h e ii is the length of the shell burning phase. For 
Si shell burning, TV ~ 400. This relatively small number 
demonstrates that the stochastic nature of the wave emis¬ 
sion process can be important during late burning phases 
(in contrast to, e.g., He burning when TV ~ 2 x 10'). Af¬ 
ter TV steps of the random walk, the core AM in each 
direction has a Gaussian distribution centered around 
zero and with standard deviation 



The magnitude of the total AM vector then has a 


Maxwellian distribution, with standard deviation <jj. 
The corresponding spin frequency also has a Maxwellian 
distribution, given by 

/(H) - . (8) 

where ctq = and I c is the moment of inertia of the 

radiative core. A little algebra yields the corresponding 
spin period distribution, 



with <t p = 27r/(7n- 

The expected angular spin frequency corresponding to 
equation 8 is 


Hex — 


luyffThell •/u 


37r 2 I c 

with a corresponding expected spin period of 


-Pex = 8/H e 


( 10 ) 

( 11 ) 


The expected spin frequency scales as H ex oc 
uj c 1/2 L c T^ u . Thus, more energetic and long-lived con¬ 
vection yields higher expected rotation rates. The short 
duration T s heii of later burning phases largely counteracts 
their increased vigor, and we find that earlier burning 
phases are generally capable of depositing more AM. 

We consider the estimate of equation 10 to be robust 
against many uncertainties associated with IGW gen¬ 
eration. It is not sensitive to the details of the IGW 
spectrum, and only relies on the fact that the spectrum 
is peaked near frequencies of u> c . Moreover, the scaling 
above shows that the value of H ex depends weakly upon 
the convective turnover frequency u> c . The main uncer¬ 
tainty is the value of the IGW flux, so we will consider 
both the optimistic and pessimistic limits of equation 1. 

The stochastic spin-up process described above will 
only occur under certain conditions. First, the core and 
burning shell must be slowly rotating, or else the stochas¬ 
tic spin-up will have a negligible effect. Second, stochas¬ 
tic spin-up can only proceed as long as H ex < oj c . If H ex 
approaches w c , wave filtering processes as described in 
Section 4 will alter the subsequent dynamics. Most of 
our estimates below have H ex <C w c , so we believe they 
are valid estimates of minimum spin rates. For C burning 
in the optimistic wave flux estimate (E oc A4 5 / 8 ), how¬ 
ever, H ex approaches w c , so this value of H ex lies near 
the maximum rotation rate achievable through stochas¬ 
tic spin-up for our stellar model. 

Moreover, stochastic spin-up can only occur if other 
sources of AM transport (e.g., magnetic torques) oper¬ 
ate on longer time scales. This could be the case during 
late burning phases when magnetic torques become in¬ 
effective (Hegcr et al. 2005; Wheeler et al. 2014). We 
can estimate a minimum magnetic coupling time be¬ 
tween core and envelope via the Alfven wave crossing 
time tA ~ fc\/~Pc/B 1 with B the approximate magnetic 
field strength. Typical NS field strengths of a few times 
10 12 G imply field strengths of ~ 10 8 G in the iron core, 
which yields ~ 5 x 10 4 s, much longer than the Si shell 
burning time (see Table 1). Although magnetic torques 
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TABLE 3 


Burning Phase 

F ex (s) 

Pex,Fe ( s ) 

Fex.NS (s) 

Pmin (s) 

Pmin,Fe (s) 

Lnin.NS (s) 

He Core Burn 

- 

- 

- 

2 x 10 5 

30 

3 X 10 -3 

C Shell Burn 

2 x 10 4 

4 X 10 2 

4 x 10 -2 

2 x 10 3 

40 

4 x 10 -3 

O Shell Burn 

2 x 10 4 

10 3 

io- 1 

- 

- 

- 

Si Shell Burn 

4 x 10 3 

2 X 10 3 

2 x 10 1 

- 

- 

- 


Approximate maximum and minimum rotation periods enforced by IGW after the burning phases shown in Figure 1 for our 
12 Mg model. The maximum rotation periods P ex are calculated from equation 11 using a conservative IGW flux, while the 
minimum rotation periods P m i n are calculated in Section 4. P ex is most plausibly set by Si shell burning because magnetic 
torques cannot compete with IGW during this phase. For each phase of evolution, we list the rotation period immediately 
after the burning phase, the corresponding rotation rate of the iron core assuming AM conservation (Fe subscript), and the 
corresponding rotation rate of a NS (NS subscript). 


TABLE 4 


Burning Phase 

F ex (s) 

Pex,Fe (s) 

F ex ,NS ( s ) 

C Shell Burn 

3 x 10 3 

50 

3 x 10 -3 

O Shell Burn 

2 x 10 3 

150 

10 -2 

Si Shell Burn 

2 x 10 3 

400 

3 x 10" 2 


Same as Table 3, but listing the values of P ex for our 35 Mg model. 
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Fig. 3. — Normalized distribution of maximum spin periods in the pre-collapse iron core, -P m ax,Fe 3 due to stochastic spin-up via IGW in 
our 12 Mq model. The left panel shows the distribution using the pessimistic wave flux (left hand side of equation 1), while the right panel 
shows the distribution using the optimistic wave flux (right hand side of equation 1). The top axis shows the corresponding maximum NS 
rotation period, -P m ax,NS 3 assuming conservation of AM during the supernova. P m ax,NS is most plausibly set during Si burning (see text), 
which corresponds to the shaded red area. Stochastic IGW spin-up leads to values of P max ns °f tens to hundreds of milliseconds, similar 
to the initial rotation rates of most pulsars. 
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may suppress stochastic spin-up during He/C/O burn¬ 
ing phases, we expect them to have a negligible impact 
during Si burning. 1 

As discussed above, the IGW may generate significant 
amounts of shear within the iron core, i.e., the core will 
be differentially rotating at CC. After collapse, we expect 
the AM to be redistributed within the subsequent NS 
on a short timescale, e.g., on a NS Alfven time scale 
Ot ~ 1 minute. Although differential rotation may persist 
through CC and the early explosion phase, we expect 
nearly rigid rotation by the time the NS spin can be 
observed many years later. 

Table 3 lists maximum expected rotation periods P ex 
calculated from equation 11 after late shell burning 
phases. We evaluate the core moment of inertia I c at 
a mass coordinate of 1.3 Mq for these estimates. We 
also list the corresponding rotation period P e x,Fe f° r the 
pre-collapse iron core, calculated from the value of I c just 
before collapse when the mass coordinate at 1.3M© has a 
radius of 1500 km. Additionally, we list the correspond¬ 
ing spin period P ex ,NS for a NS with Mns = 1.3Mq, 

Pns = 12 km, and /ns = 0.25Mns/?ns- 

Figure 3 shows the distribution in maximum spin pe¬ 
riod of the pre-collapse iron core, P ma x,Fe- We have plot¬ 
ted the values of P ma x,Fe, if the core spin rate is set during 
C, O, or Si burning. We find that C, O, and Si burn¬ 
ing generate maximum iron core rotation periods in the 
range of tens to thousands of seconds, depending on the 
burning phase and efficiency of IGW generation. Si burn¬ 
ing most plausibly sets the pre-SN conditions because it 
is the last convective burning phase before CC and is 
least likely to be affected by magnetic torques. For the 
conservative wave flux, we find stochastic spin-up during 
Si burning leads to 400 s < P max , Fe ^5x 10 3 s. The cor¬ 
responding NS rotation rate is 40 ms < Pns ^5 500 ms. 
Hence, we find that very slow core rotation rates, as sug¬ 
gested by Spruit & Phinncy (1998), are unlikely. Nor do 
we expect that there is a population of NSs born with 
very long spin periods, P > 2 s, at least from progenitors 
with ZAMS mass 10 Mq < M < 20Mq. 

The distribution of NS spin periods shown in Figure 
3 appears broadly consistent with those inferred for 
young NSs (Faucher-Giguere & Kaspi 2006; Popov et al. 

1 It is possible that the stochastic spin-up process saturates due 
to back-reaction on the convective shell, which gains AM opposite 
to the AM deposited by IGW in the core. The induced spin of the 
convective shell may alter wave generation such that the AM of sub¬ 
sequent wave packets is not randomly oriented and the AM of the 
core does not undergo a purely random walk. This effect depends 
on uncertain details such as the rotation profile at the core-shell 
interface and the effect of the Coriolis force on wave generation. 
However, we may guess that in the limit of Q s | ie ii <C ui c , the AM of 
each wave packet could obtain a non-stochastic component of order 
m\ lilsheii/tUc. In this case, the stochastic build-up of AM within 
the core would saturate when the net deposition of stochastic AM 
is comparable to the non-stochastic AM deposition, which occurs 
when VN ~ Nm£l s h e \\/uj c . The expected rotation rate of the core 

after stochastic spin saturation is Sl e x ~ \J 2?r/shell®/UI C / Ic- We 
find that stochastic spin-up may saturate during C/O burning, 
leading to P ex p e ~ 2 x 10 3 s for each of these stages. Interest¬ 
ingly, this spin period is similar to the value of P ex ,Fe generated 
by stochastic spin-up during Si burning. However, the saturation 
state is unlikely to be reached during Si shell burning, and so our 
calculations for Si burning remain unchanged. Therefore, we find 
the value of P ex set by Si burning to be reasonable, regardless of 
the details of stochastic spin saturation. 


2010; Gullon et al. 2014). We are therefore tempted to 
speculate that the stochastic spin-up scenario described 
above may be the dominant process setting the spin 
rates of newly born NSs. If so, this scenario predicts that 
the rotation rate and direction of the NS is uncorrelated 
with the rotation of the envelope of the progenitor 
star, in contrast to any sort of magnetic spindown 
mechanism. However, there are several caveats to keep 
in mind. First, the scenario presented above can only 
proceed if the core is initially very slowly rotating, which 
requires efficient magnetic/IGW core spin-down (see 
Section 4) to occur before Si burning. Second, the NS 
rotation rate may be changed during the supernova, by 
fallback effects, or by the r-rnode instability (Andersson 
1998; Andersson ct al. 1999, see Ott 2009 for a review). 
Finally, there is a considerable amount of uncertainty in 
the IGW energy flux. Since the minimum core rotation 
rate set by IGW is proportional to the wave energy flux 
(which is uncertain at an order of magnitude level), 
there is an equal amount of uncertainty in the induced 
rotation rates. 

To understand how stochastic IGW spin-up proceeds 
in different types of massive stars, we have performed the 
procedure above for a 35 Mq ZAMS model (described 
in Appendix A), whose late shell burning properties are 
listed in Table 2. The expected spin periods P ex (equa¬ 
tion 11) are listed in Table 4, and the distribution of 
expected periods is plotted in Figure 4. In general, the 
spin periods of our 35 M 0 model are smaller by a fac¬ 
tor of a few compared to our 12 Mq model. The reason 
is that the later burning phases (the O burning phase 
in particular) are more vigorous (higher L c ) and have a 
shorter duration (lower T s h e ii) in higher mass stars with 

larger He cores. Since H ex oc T sh ' eU I/ c , the increased con¬ 
vective luminosity wins out, generating higher spin rates 
in larger mass stars. It is therefore possible that higher 
mass stars give birth to more rapidly rotating NSs (unless 
they form black holes instead). For the optimistic wave 
flux (right panel) or for stochastic spin-up via O shell 
burning, NS spin periods on the order of milliseconds 
could be generated via stochastic IGW spin-up. We note 
that the shorter duration O burning phase of higher mass 
stars makes it more likely that stochastic spin-up during 
O shell burning can be preserved until CC. However, as 
the O shell burning lifetime /’shell ~ 10 days may still be 
longer than an Alfven crossing time, magnetic torques 
may suppress stochastic spin-up during C/O burning as 
discussed above. 

Unfortunately, we cannot hope to map out the pre¬ 
collapse spin rates of all massive stars in this paper. Dif¬ 
ferent progenitors may yield qualitatively different re¬ 
sults. For instance, electron-capture SNe progenitors, 
whose final stages burning stages are quite different from 
the CC case, may have different core spin rates. Very 
massive stars (M > 50Mq) and stars altered by binary 
evolution will also make interesting targets for future 
studies. 

4. CORE SPIN-DOWN BY INTERNAL GRAVITY WAVES 

In Section 3, we discussed the spin-up of a non-rotating 
core during the final shell burning phases before CC. 
However, stars are born rotating, and in the absence of 
AM transport the core will spin-up as it contracts during 
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Fig. 4. — Same as Figure 3, but for our 35 Mq progenitor. Larger burning luminosities in this model lead to larger stochastic spin-up 
rates and smaller values of P max . O burning may be more important for this model, but Si burning most likely sets the value of P max (see 
text). 


stellar evolution. Our goal here is to determine whether 
IGW generated during shell burning phases can prevent 
the spin-up of the core due to contraction. Because we 
are interested in understanding the effects of IGW on a 
relatively rapidly rotating stellar core, the random walk 
arguments of the previous section do not apply. Instead, 
as we now explain, the key physics is the preferential 
absorption of certain waves due to differential rotation- 
induced wave damping. 

Over long time scales, we expect turbulent convec¬ 
tion to generate prograde and retrograde waves in nearly 
equal quantities, so that the net AM flux imparted to 
IGW is nearly zero. However, differential rotation can set 
up powerful wave filtration mechanisms (see F14), which 
filter out either prograde or retrograde waves. Consider 
waves of angular frequency ui and azimuthal number m 
that propagate across a radial region of thickness Ar, 
whose endpoints have angular spin frequencies that dif¬ 
fer by an amount AH. If AH > u/m, the waves will 
encounter a critical layer within the region, and will be 
absorbed. Only waves of opposing AM will penetrate 
through the layer; therefore, rapidly rotating regions will 
only permit influxes of negative AM, which will act to 
slow the rotation of the underlying layers. 

IGW can therefore limit differential rotation to a max¬ 
imum amplitude AH max ~ w/m, provided the IGW AM 
flux is large enough to change the spin rate on time 
scales shorter than relevant stellar evolution times. In 
the case of a rapidly rotating core (which has contracted 
and spun-up) surrounded by a slowly rotating burning 
shell, we may expect a maximum core rotation rate of 
~w, provided waves of this frequency can propagate into 
the core. This maximum rotation rate assumes |m| = 1 
waves dominate the AM flux, the actual rotation rate 


could be smaller if \m\ > 1 waves have a substantial im¬ 
pact. Thus we define a maximum core rotation rate 

flmax ^ ; (12) 

where is the characteristic frequency of waves that 
dominate AM transport and are able to penetrate into 
the core. 

The wave frequency w* which dominates AM trans¬ 
port is determined by the IGW frequencies generated 
by a convective shell, and by the subsequent propaga¬ 
tion and dissipation of those waves. In the absence of 
wave damping, w* ~ w c because these waves dominate 
the energy/AM flux. However, radiative diffusion pref¬ 
erentially damps low frequency IGW because they have 
shorter radial wavelengths and slower group velocity. As 
low frequency IGW damp out, the value of w* shifts to 
larger frequencies at larger depths below the convective 
zone. Although radiative dissipation of ingoing IGW 
is negligible from C shell burning onward, we find it is 
important for IGW generated by the surface convection 
zone during the core He burning phase. 

To calculate the appropriate value of w*, we use the 
same methods as F14. Upon generation, the IGW carry 
an energy flux Eq and AM flux Jo, and have a frequency 
spectrum which is initially peaked around u c . We assume 
the wave spectrum has a power law fall off at higher 
frequencies such that 

(13) 

dui uj c \uj c J 

where a is the slope of the frequency spectrum, which is 
somewhat uncertain. As in F14, we expect a spectrum 
slope in the range 3 < a < 7, with a fiducial value of 
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a = 4.5 as found by Kumar et al. 1999; Talon et al. 2002. 

F14 show that radiative damping leads to a value of 
to* of 


to* (r) = max 


dr 


X 3 / 2 N^NK \ 1/4 

r 3 J 


(14) 


The corresponding AM flux carried by waves of to ~ to* is 


j*(r) 


to* (r) 
to c 


—a 


Jo- 


(15) 


In equation 14, r c is the radius of the inner edge of the 
convective zone, A = 1(1 + 1), l is the angular index of 
the wave (which corresponds to its spherical harmonic 
dependence, Yj m ), Nt is the thermal part of the Brunt- 
Vaisala frequency, and K is the thermal diffusivity. In 
what follows, we focus on / = 1 waves because they have 
the longest damping lengths and will dominate the AM 
flux when the waves are heavily damped. Moreover, fo¬ 
cusing on l = 1 waves allows us to estimate maximum 
spin frequencies, although slower spin frequencies can be 
obtained when higher values of l and m contribute to the 
AM flux. 

The top panel of Figure 5 shows the value of to*(r) dur¬ 
ing different burning phases for our 12 M 0 model. For 
C, O, and Si burning, radiative diffusion is negligible and 
to* ~ to c everywhere. Thus, our calculations are insensi¬ 
tive to uncertainties in the wave spectrum during these 
stages, to* increases near the center of these models due 
to non-linear breaking (see Appendix B). During He core 
burning, however, radiative diffusion is very important, 
causing the value of to* to increase by a factor of ~ 10 as 
IGW propagate inward. This decreases the AM carried 
by the waves and the minimum rotation period they can 
enforce. 

The AM deposited by IGW is only significant if it is 
larger than the amount of AM contained within the core 
of the star. A typical massive star has a zero-age main se¬ 
quence equatorial rotation velocity of v ro t ~ 150 km s ^ 1 
(de Mink et al. 2013), corresponding to a rotation period 
of -Pms ~ 1.5 d for our stellar model. Using this rotation 
rate, we calculate the AM Jo (M) contained within the 
mass coordinate M(r), given rigid rotation on the main 
sequence. In the absence of AM transport, this AM is 
conserved, causing the core to spin up as it contracts. 
Of course, magnetic torques may extract much of this 
AM, so Jo represents an upper limit to the AM contained 
within the mass coordinate M ( r ). Both Jo and the corre¬ 
sponding evolving rotation profiles are shown in Figure 5. 
We also plot the approximate AM Jns contained within 
a NS rotating at Pns = 10 ms, which is more than two 
orders of magnitude smaller than the value of Jo within 
the inner 1.4M 0 . 

The AM capable of being removed from below a radius 
r via IGW launched during a burning phase is 


Jex — (r ) Jshell ■ ( 16 ) 

Here, r) is an efficiency factor that accounts for the fact 
that only some of the wave flux is in low l retrograde 
waves capable of depositing negative AM in the core. We 
estimate the value of 77 as follows. For waves with to ~ to c 
launched from a thick convective zone with r ~ H, the 
energy spectrum described in Kumar et al. (1999); Talon 



M(M 0 ) 


Fig. 5.— Top: Angular frequency a;* of waves that dominate 
AM transport within our 12 M© model at different phases of evolu¬ 
tion. The dashed vertical lines mark the bottom of the convective 
zone from which the IGW are launched, while the horizontal ex¬ 
tent of each line marks the cavity in which the IGW propagate. 
We have truncated the curves for He burning at the location of 
the convective core, into which the IGW cannot propagate. The 
value of cu* increases inward for the He burning model because 
lower frequency waves are damped by radiative diffusion, although 
this effect is negligible for C burning and beyond. Middle: AM 
Jo (thick black line) contained within the mass coordinate M(r) 
of our model rotating with a period of Po ~ 1.5 d on the ZAMS, 
while the orange line shows the approximate AM Jns contained 
within M(r) for a NS rotating at Pns = 10 ms. The shaded re¬ 
gions indicate the AM J ex (equation 16) that can be extracted by 
IGW during each burning phase. The shaded regions are bounded 
by lines calculated with a pessimistic and optimistic estimate for 
IGW fluxes (see equation 1). The hatched region for He burning 
is calculated using the pessimistic IGW flux for a steep IGW fre¬ 
quency spectrum (a = 5.5, lower bound) and a shallow spectrum 
(a = 3.5, upper bound). IGW can significantly slow the spin rate 
of the progenitor in regions where J ex > Jo- Bottom: Spin peri¬ 
ods Po of our model in the absence of AM transport. As the star 
evolves, the value of Po decreases in the contracting core. We have 
also plotted the approximate minimum spin periods P m i n which 
may be enforced by IGW during the core He burning and C-shell 
burning phases. 


ct al. (2002); Talon & Charbonnel (2005) has an approx¬ 
imate spectrum dE/dl oc le~ l . With this spectrum, the 
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energy emitted in l = 1 waves is approximately one third 
of the total energy flux. These waves can have azimuthal 
numbers m = —1, m = 0, or m = 1, therefore we expect 
approximately one third of l = 1 waves are retrograde 
waves with m = — 1. Hence, we find ?7 = 0.1 to be a 
reasonable estimate for the energy flux emitted in low 
degree retrograde waves capable of propagating into the 
core. However, we caution that there is a great deal of 
uncertainty in the wave spectrum, and these numbers 
should be viewed only as order of magnitude estimates. 

Figure 5 shows both a pessimistic and optimistic es¬ 
timate for J ex , corresponding to the left and right-hand 
sides of equation 1, respectively. We find that the values 
of Jex are comparable to Jo for waves emitted during He 
core burning and C shell burning. This implies that IGW 
emitted during these phases may be able to significantly 
spin down the cores of massive stars. During core He 
burning, the inner ~2Mq is convective and IGW cannot 
propagate into it, hence, it will only be spun down if it 
is coupled (via its own IGW or via magnetic torques) to 
the radiative region above it. We have also considered 
a range of IGW spectra corresponding to 3.5 < a < 5.5 
in equation 13. For steep spectra (a = 5.5), there is not 
enough power in high frequency IGW to allow them to 
spin down the core. For shallow spectra (a = 3.5), IGW 
likely can spin down the core. More sophisticated anal¬ 
yses (predicated on a better understanding of the IGW 
spectrum) are required for a robust conclusion. 

During O and Si shell burning, we find that IGW most 
likely cannot remove the AM contained within the core, 
if the core retains its full AM from birth. This does not 
imply IGW have no effect, as the value of J ex for O/Si 
burning is larger than Jns (the typical AM content of 
a fairly rapidly rotating NS). Therefore, if the core has 
been spun down by IGW or magnetic torques during pre¬ 
vious burning phases, IGW during late burning phases 
may be critical in modifying the core spin rate (as dis¬ 
cussed in Section 3). 

If IGW are able to spin down the core during He core 
burning or C shell burning, this entails a minimum pos¬ 
sible core rotation period P m i n = 27r/w*(r) at the end 
of these phases. The bottom panel of Figure 5 plots 
the value of P m ; n , in addition to the rotation profile Pq 
corresponding to the AM profile Jo that would occur in 
the absence of AM transport. If IGW are able to spin 
down the cores, the minimum rotation periods are 10- 
100 times larger than those that would exist without AM 
transport. Thus, IGW may significantly spin down the 
cores of massive stars. We also note that the value of 
w* is not sensitive the the value of a, so P m in is insensi¬ 
tive to the IGW spectrum, as long as IGW have enough 
power to spin down the core. Table 1 lists the values of 
Pnin corresponding to He and C burning, as well as corre¬ 
sponding minimium spin periods for the pre-collapse iron 
core (P m in,Fe) and for the neutron star remnant (P m in,Ns) 
given no subsequent AM transport. The minimum NS 
rotation period P m i nj NS we calculate is on the order of 
milliseconds, which is shorter than that inferred for most 
newly born NSs. Therefore either IGW spin-down is sig¬ 
nificantly more effective than our estimates, or (perhaps 
more likely) magnetic torques are responsible for spin¬ 
ning down the cores of massive stars. 


5. DISCUSSION AND CONCLUSIONS 

We have demonstrated that convectively generated in¬ 
ternal gravity waves (IGW) in massive stars are capa¬ 
ble of redistributing angular momentum (AM) on short 
time scales. We have focused primarily on the effects of 
IGW generated during late stages of massive star evolu¬ 
tion (He burning and later) for typical neutron star (NS) 
progenitors (1OM 0 < M < 20M 0 ). It may seem sur¬ 
prising that AM transport via IGW can act on the short 
stellar evolution timescales of massive stars nearing core 
collapse (CC). However, the huge convective luminosi¬ 
ties inside evolved massive stars ensure large fluxes of 
IGW (QS12, SQ14) that can transport energy and AM 
on short timescales. We therefore encourage efforts to 
incorporate the effects of IGW in stellar evolution codes 
focusing on the final stages of massive star evolution. 

During He/C burning, inwardly propagating IGW 
launched from convective shells may slow down the core 
to much slower spin rates than would be obtained in the 
absence of other AM transport mechanisms, although 
the result depends on the IGW spectrum and a robust 
conclusion remains elusive. If IGW do have an impact, 
spin-down during core He burning may slow the outer ra¬ 
diative core to minimum spin periods of P m i n ,He ~ 2 days 
in our 12 M 0 model. The inner He burning core will also 
be spun down if it is strongly coupled with the outer core 
via IGW or magnetic torques. IGW launched during C 
shell burning may also be able to substantially slow the 
spin of the core. These convective phases plausibly lead 
to pre-collapse iron cores with minimum rotation periods 
Pmin.Fe > 30 s, corresponding to initial NS rotation peri¬ 
ods of P m in,NS ^ 3 ms. The rotation periods listed above 
are minimum periods for our stellar model. Calculations 
of rotation rates including magnetic torques (Heger et al. 
2005; Wheeler et al. 2014) typically yield rotation periods 
several times larger. Magnetic torques may therefore be 
the dominant AM transport mechanism responsible for 
extracting AM from massive stellar cores, although it is 
possible that both mechanisms play a significant role. 

Stochastic influxes of IGW during late burning phases 
can also lead to the spin-up of an otherwise very slowly 
rotating core. This occurs in the case of very efficient 
prior core spin-down via IGW/magnetic torques. Such 
efficient core spin-down is not unreasonable, especially 
given that the cores of low mass red giant stars rotate 
slower than can be accounted for using existing prescrip¬ 
tions for hydrodynamic mechanisms or magnetic torques 
via the Tayler-Spruit dynamo (Canticllo et al. 2014). It is 
thus quite plausible that massive star cores are efficiently 
spun down via IGW/magnetic torques, after which they 
are stochastically spun up via IGW launched during O/Si 
burning. This process is similar to previous theories 
of IGW-induced spin alteration of massive stellar atmo¬ 
spheres (Rogers et al. 2012), and stochastic spin-up of 
proto-NSs during supernovae (Spruit & Phinney 1998). 
The interesting feature of IGW spin-up during the final 
stages of massive star evolution is that it may occur on 
time scales shorter than a core Alfven crossing time and 
therefore can operate without suppression from magnetic 
torques. 

If this mechanism determines the core spin rate before 
death, it predicts a Maxwellian distribution in spin fre¬ 
quency, with typical iron core spin periods of 400 s < 
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Ppe 1$ 5 x 10 3 s. We thus find it unlikely that magnetic 
torques can enforce very large pre-collapse spin periods 
as claimed by Spruit & Phinney (1998). Additionally, 
we speculate that the stochastic spin-up process is rel¬ 
atively insensitive to binary interactions or winds that 
have stripped the stellar envelope, as long as these pro¬ 
cesses do not strongly modify the core structure and late 
burning phases. We also express a word of caution, as Si 
burning is notoriously difficult for stellar evolution codes 
to handle, and the properties of Si burning produced by 
our MESA evolutions have large associated uncertainties. 
The rough energy and AM fluxes in convectively excited 
waves are, however, reasonable at the order of magnitude 
level. 

If AM is conserved during the supernova, stochastic 
IGW spin-up entails NS birth periods of 40 ms < Pns < 
500 ms, albeit with significant uncertainty. These es¬ 
timates are comparable to initial spin periods of some 
young NSs (Lai 1996; Gotthelf et al. 2013), and to the 
broad inferred birth spin period distribution of Pns 
500 ms for ordinary pulsars (Fauchcr-Gigucre & Kaspi 
2006; Popov et al. 2010; Gullon et al. 2014). Therefore, 
stochastic IGW spin-up could be the dominant mecha¬ 
nism in determining the rotation periods of pre-collapse 
SN cores and newborn NSs. In this scenario, there is 
little or no correlation between the spin of the progeni¬ 
tor and the spin of the NS it spawns. Although torques 
during the supernova may modify the spin rate of the 
NS, they would have to be very finely tuned to erase the 
stochastic spin-up occurring during shell burning. Any 
sort of purely frictional spin-down processes would likely 
slow the NS to rotation periods larger than typically in¬ 
ferred for young NSs. 

We have also investigated stochastic IGW core spin- 
up for a 35 M 0 model. We find that the increased late 
time burning luminosities of more massive stars leads to 
more effective stochastic spin-up, creating spin rates a 
factor of a few higher than those listed above. For opti¬ 
mistic wave fluxes, the initial rotation rates of NSs born 
from massive progenitors may be as short as several mil¬ 
liseconds. Even larger rotation rates could be possible if 
magnetic torques are weak enough to allow the stochastic 
core spin-up generated during O shell burning to persist 
until CC. 

We have quantitatively considered the implications of 
IGW AM transport in two distinct limits. First, we ne¬ 
glect AM transport by magnetic fields and consider the 
limit in which the core is rotating much faster than the 
convective shell because of AM conservation during core 
contraction. In this case, IGW emitted from convective 
shells propagate into the radiative core and may be able 
to substantially slow its rotation, enforcing a maximum 
rotation rate. The second limit we consider is when the 
stellar core has been efficiently spun-down via magnetic 
coupling to the envelope and/or IGW in earlier phases 
of stellar evolution. In this case, we have shown that 
the stochastic influx of AM via IGW from shell burning 
leads to a spin-up of the stellar core and a minimum core 
rotation rate. Taken together, these limits enforce iron 


core rotation periods 30 s < Pp e < 5 x 10 3 s and initial 
NS rotation periods of 3 ms < Pns 500 ms. We expect 
these limits to be robust against many uncertain factors 
in massive star evolution, e.g., birth spin rate, mass loss, 
mixing, and the effects of magnetic fields. 

There is ample evidence that some CC events occur 
with rapidly rotating cores. In particular, long GRBs al¬ 
most certainly require a rapidly rotating central engine 
(Woosley 1993; Yoon et al. 2006; Woosley & Heger 2006; 
Metzger et al. 2011), and the picture advanced above 
must break down in certain (although somewhat rare) 
circumstances. It is not immediately clear what factors 
contribute to the high spin rate in GRB progenitors, as 
our analysis was mostly restricted to “typical” effectively 
single NS progenitors with 1OM 0 < M < 20M 0 , which 
explode to produce type-IIp supernovae during a red su¬ 
pergiant phase (see e.g. Srnartt 2009). We speculate that 
GRB progenitors (if occurring in effectively single star 
systems) have never undergone a red supergiant phase, 
as torques via magnetic fields and/or IGW are likely to 
spin down the helium core by coupling it with the huge 
AM reservoir contained in the slowly rotating convective 
envelope. Our preliminary examination of more mas¬ 
sive models indicates that stochastic spin-up may lead to 
larger spin rates, but is unlikely to generate very rapid 
rotation throughout the entire He core. A third possi¬ 
bility is that spin-up via mass transfer/tidal torques in 
binary systems is required for GRB production (( lanticllo 
et al. 2007). A merger or common envelope event after 
the main sequence could also remove the extended con¬ 
vective envelope and prevent it from spinning down the 
core. 

The population of massive stars approaching death 
is complex, and factors such as initial mass, rotation, 
metallicity, binarity, magnetic fields, overshoot, mixing, 
winds, etc., will all contribute to the anatomy of ag¬ 
ing massive stars. We have argued that AM transport 
via convectively driven IGW is likely to be an impor¬ 
tant factor in most massive stars. But it is not imme¬ 
diately obvious how this picture will change in different 
scenarios, e.g., electron capture supernovae, very massive 
(M > 50M 0 ) stars, interacting binaries, etc. We hope 
to explore these issues in subsequent works. 
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APPENDIX 
STELLAR MODELS 

Our stellar models are constructed using the MESA stellar evolution code (Paxton et al. 2011, 2013), version 6794. 
The models are run using the following inlist controls file. 


\& star_job 

kappa_file_prefix = ’gs98’ 
/ ! end of star_job namelist 


\& controls 

initial_mass = 12. 
initial_z = 0.020 

sig_min_factor_for_high_Tcenter =0.01 
Tcenter_min_for_sig_min_factor_full_on = 3.2d9 
Tcenter_max_for_sig_min_factor_full_off = 2.8d9 
logT_max_for_standard_mesh_delta_coeff = 9.0 
logT_min_for_highT_mesh_delta_coeff = 10 
delta_Ye_highT_limit = ld-3 
okay_to_reduce_gradT_excess = .true. 
allow_thermohaline_mixing = .true. 
thermo_haline_coeff = 2.0 
overshoot_f_above_nonburn = 0.035 
overshoot_f_below_nonburn = 0.01 
overshoot_f_above_burn_h = 0.035 
overshoot_f_below_burn_h = 0.0035 
overshoot_f_above_burn_he = 0.035 
overshoot_f_below_burn_he = 0.0035 
overshoot_f_above_burn_z = 0.035 
overshoot_f_below_burn_z = 0.0035 
RGB_wind_scheme = ’Dutch’ 

AGB_wind_scheme = ’Dutch’ 

RGB_to_AGB_wind_switch = ld-4 
Dutch_wind_eta = 0.8 
include_dmu_dt_in_eps_grav = .true. 
use_Type2_opacities = .true. 
newton_itermin = 2 
mixing_length_alpha = 1.5 
MLT_option = ’Henyey’ 
allow_semiconvective_mixing = .true. 
alpha_semiconvection = 0.01 
mesh_delta_coeff = 1. 
varcontrol_target = 5d-4 
max_allowed_nz = 10000 
mesh_dlog_pp_dlogP_extra = 0.4 
mesh_dlog_cno_dlogP_extra = 0.4 
mesh_dlog_burn_n_dlogP_extra = 0.4 
mesh_dlog_3alf_dlogP_extra = 0.4 
mesh_dlog_burn_c_dlogP_extra =0.10 
mesh_dlog_cc_dlogP_extra =0.10 
mesh_dlog_co_dlogP_extra =0.10 
mesh_dlog_oo_dlogP_extra =0.10 
velocity_logT_lower_bound=9 
dX_nuc_drop_limit=5d-3 
dX_nuc_drop_limit_at_high_T = 5d-3 ! 
screening_mode = ’extended’ 
max_iter_for_resid_toll = 3 
tol_residual_norml = ld-5 
tol_max_residuall = ld-2 
max_iter_for_resid_tol2 = 12 
tol_residual_norm2 = ld99 
tol_max_residual2 = ld99 
min_timestep_limit = ld-12 ! (seconds) 
delta_lgL_He_limit =0.1 ! 
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dX_mic_drop_max_A_limit = 52 
dX_niic_drop_min_X_limit = ld-4 
dX_nuc_drop_hard_limit = ld99 
delta_lgTeff_limit = 0.5 
delta_lgL_limit = 0.5 
delta_lgRho_cntr_limit = 0.02 
T_mix_limit = 0 
/ ! end of controls namelist 


The most important feature of this model is that it contains significant convective overshoot, especially above 
convective zones. It is non-rotating, thus there is no rotational mixing. 

Just before core O burning, we change to a 201-isotope reaction network: 

change_net = .true. 
new_net_name = , mesa_201.net’ 

Although our choices affect details of the model (e.g., He core mass), the general features of our model are robust. 
It always explodes as a red supergiant. It always undergoes convective core C burning, followed by shell C burning, 
core O/Ne burning, shell O burning, core Si burning, shell Si burning, and then CC. The approximate convective 
properties (as described by MLT) are not strongly affected by model parameters. Since these properties are most 
important for IGW AM transport, we argue that the general features of IGWs described in this work are fairly robust 
against uncertain parameters in our massive star models. 

Differences in the inlist for our 35 Mq model are 


&star_job 

relax_initial_Z = .true. ! gradually change abundances, reconverging at each step. 
new_Z = 2d-3 

/ ! end of star_job namelist 

&controls 

initial_mass = 35. 

okay_to_reduce_gradT_excess = .true. IMLT++ 

overshoot_f_above_nonburn = 0.015 
overshoot_f_below_nonburn = 0.005 

overshoot_f_above_burn_h = 0.015 
overshoot_f_below_burn_h = 0.005 

overshoot_f_above_burn_he = 0.015 
overshoot_f_below_burn_he = 0.005 

overshoot_f_above_burn_z = 0.015 
overshoot_f_below_burn_z = 0.005 


varcontrol_target = 8d-4 
/ ! end of controls namelist 


Notable features of this model are its lower metallicity and the use of MLT++ Paxton et al. (2013) to evolve the 
star through super-Eddington phases of evolution. 

NON-LINEAR DAMPING 

IGW will overturn and break, leading to local energy/AM deposition, if they obtain sufficiently non-linear amplitudes. 
Here we estimate those amplitudes and the AM flux that can be carried toward the center of the star as waves are 
non-linearly attenuated. For traveling waves in the WKB limit, it is well known that the radial wave number is 

VXN 

ruj 


k 


(Bl) 
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It is straightforward to show that the radial displacement £ r associated with IGW of frequency <x> carrying an energy 
flux E is 

VXE 4 1/2 


|£r| = 

The waves become non-linear and break when 

I kr^r | 


■ Y ." 

X 3 / 2 NE ] 1/2 


pr 5 w 4 


1. 


(B2) 


(B3) 


In the absence of damping, FI is a conserved quantity. Therefore, waves become more non-linear as they propagate 
into regions with larger N, lower density, or smaller radius. In our problem, the geometrical focusing (i.e., the r- 
depdendence) is the most important feature of equation B3, and causes waves to non-linearly break as they propagate 
inward. Note also the oj~ 2 dependence of equation B3, which causes low frequency waves to preferentially damp. 

Equation B3 entails there is a maximum energy flux that can be carried by waves of frequency w, 


E max 


A 2 pr 


\3/ 2 N 


(B4) 


for waves that non-linearly break when |fc r £ r | = A ~ 1. We find that IGW of frequency w ~ w c become strongly 
non-linear as they propagate inward during C/O/Si shell burning. Therefore, most of the IGW energy/AM will be 
deposited within the core. Where the waves are highly non-linear, the waves which dominate the energy flux are those 
which are on the verge on breaking. To determine their frequency, we use the frequency spectrum of equation 13 to 
find 


Er, 


Solving equation B5 yields the wave frequency which dominates energy transport, 


EA — 


w* 
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-l/(a+3) 
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(B6) 


We expect frequency spectra with slopes somewhere near 3 < a < 7. Therefore the exponent in equation B6 is quite 
small, and in most cases, w* does not increase to values much larger than w c . 

Substituting equation B6 back into equation B4 allows us to solve for the energy and AM flux as a function of radius 
due to non-linear attenuation. The result is 


J * 


w*(r) 


Jo 


Jo 


( A 2 pr 5 u> 4 
\X 3 / 2 NE q 


a/(o+3) 


Jo 


(B7) 


During C/O shell burning, we find that radiative diffusion only slightly damps waves near the shell burning convective 
zone, while non-linear breaking damps waves near the center of the star. In this case, we first calculate w* and its 
corresponding energy flux E, via equations 13 and 14. We then substitute the value of E 1 * for E 0 in equation B6. The 


«(Q 


Jo 


appropriate value of w* is then w* = max[Eqn. 14, Eqn.BG] . The corresponding AM flux is J* ~ 

The cores of massive stars nearing death cool primarily through neutrino emission, so it is not unreasonable to think 
that waves may be damped via neutrino emission. We calculate neutrino energy loss rates in the same manner as 
Murphy ct al. (2004). We find that neutrino damping time scales are always longer than the wave crossing timescale 


dr 

v a 


(B8) 


where the IGW radial group velocity is v g = rui 2 /(\/XN). This is not surprising, as Murphy et al (2004) found 
neutrino growth/damping rates were slower than stellar evolution time scales. Neutrino damping can therefore be 
safely ignored. 


WAVE-FLOW INTERACTION 


It is well known that the interaction between flows (e.g., differential rotation) and IGW can strongly affect IGW 
propagation and dissipation. Here we show that these interactions will not prevent the stochastic spin-up process 
(Section 3) from occurring. 

One reason that wave-flow interactions will not prevent stochastic spin-up is that the waves are very weakly damped 
by radiative diffusion during late burning stages. As shown in Section 4, radiative damping times are long. Nonetheless, 
the background flows are unstable in the presence of IGW in the sense that IGW will amplify very small amounts of 
shear. From equation 13 of F14, the timescale on which waves amplify shear due to radiative diffusion is 
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pr 4 u;*L d 

J * 


(Cl) 






















Waves in massive stars 


15 



M(M 0 ) 


Fig. 6. — Growth time for the development of shear, t gTOW , due to radiative diffusion in the radiative cores of our 12 M@ model. For 
O/Si burning, t gro w T’sheii? and it is unlikely that radiative diffusion will cause the IGW to generate significant amounts of shear. 


where Ld is the wave damping length due to radiative diffusion. In the Sun, this timescale is quite short. However, 
during late burning stages in the cores of massive stars, it typically exceeds the remaining lifetime of the star. We plot 
fg r ow below in Figure 6. During 0/Si burning, f grow is much longer than the life-time of the star. During C-burning, 
it is comparable to the remaining lifetime, and therefore IGW may create large amounts of shear during C-burning. 

Because i grow is long, IGW will not be able to amplify small amounts of shear into critical layers that absorb 
subsequent waves and prevent IGW propagation into the core. Therefore, we find that the stochastic spin-up process 
outlined in Section 3 will not be prevented by IGW-flow interaction. 
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